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ABSTRA.CT 



Waves in the surf zone are a highly nonlinear process 
which is evident by the appearance of secondary waves. The 
secondary waves appear as strong peaks in the period PDFs 
corresponding to the first harmonic of the peak of the wave 
spectrum. The strong first harmonic period peak is also 
reflected in the highly correlated height and velocity PDFs. 

Breaking waves may be viewed as a gradation between two 
extremes. The period, height and velocity PDFs for spilling 
breakers generally were unimodal, whereas for plunging 
breakers strong bimodality was found. 

The joint probability density functions for periods and 
heights of the breaking waves show high correlation ( 0 . 60 - 
0 . 80 ) which says that greater wave periods are associated 
with larger breal'.er heights. The joint PDFs of period and 
particle velocity, and velocity and height, suggest that the 
maximum onshore particle velocities are correlated with both 
the wave periods and wave heights. 

Due to the high probability of the secondary waves, the 
mean wave period for breakers is a poor descriptor of the 
average period of the offshore incident waves. 
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I . INTRODUCTION 



Waves in deep water can generally be represented to a 
first approximation as a Gaussian process having zero skew- 
ness and kurtosis equal to three. The random sea surface 
is viewed as the linear superposition of statistically inde- 
pendent spectral components which is adequately described 
by the energy density spectrum. The restriction of statis- 
tical independence implies the v/ave components act as free 
waves and there is no coupling or interaction between 
various spectral components. 

As the waves begin to shoal, the crests become more 
peaked, and the troughs become more elongated; this results 
in a positive value of the third statistical moment, i.e., 
the probability density distribution of sea surface eleva- 
tions becomes skewed to the right. The skewness is a result 
of phase coupling of the wave components, which is indica- 
tive of the increasing importance of the non-linearities 
of the system. This initially weakly non-linear system 
found in shallow water may be adequately described by a 
Gram-Charlier probability distribution (Longuet-Higgins , 
1963). 

Near breaking waves become asymmetrical about the crest 
with a continually steepening of the wave face. Miller 
(1976) has shown that a jet is then formed at the wave 
crest at breaking, which extends forward and then down 
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until it closes V7ith the forward face of the breaking wave 
form. The asymmetry about the wave crest yields a change 
in the fourth moment, or kurtosis , which is also indicative 
of the increasing non-linearities. Thornton et al. (1976) 
have noted that secondary X7ave formation becomes apparent 
on the trailing face of the wave form, due to transfer of 
energy from the primary frequency to its harmonics which 
also points to a strongly nonlinear system. 

This thesis uses data collected as described by Galvin 
(1975) and Olsen (1977) plus additional data collected by 
the author. The experiments vjere designed to measure the 
sea surface elevation and the onshore particle velocities 
by means of capacitance wave gauges and electromagnetic 
flowmeters, respectively. The experiments were conducted 
at a variety of beaches in order to include wave types of 
spilling, plunging, and collapsing forms. The characteris- 
tics of the breaking wave forms are a function of the off- 
shore wave spectrxim and the nearshore bottom slope and 
topography, where the increasing beach slope increases the 
rapidity of the breaking of the v?aves . 

Empirical distributions were obtained for the periods, 
wave heights, and maximum onshore particle velocities, for 
waves which had been defined by the zero-up - cros s method. 
Each of these distributions was compared with the Rayleigh 
probability density fimction (PDF). Joint PDFs of periods 
and wave heights, maximum velocities and wave heights, and 
periods and maximum velocities were calculated. 
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II. HISTORICAL BACKGROUND 



A. QUALITATIVE DESCRIPTION 

Thornton et al. (1976) have described a number of ob- 
served universal similarities for the various types of 
breakers occurring on different beaches. Normally there is 
a quick drawdown of water just before the breaker arrives, 
followed by a steep, nearly vertical leading edge. Miller 
(1976) suggests that subsequent to the development of the 
leading vertical edge, a jet forms at the crest. This jet 
progressively extends itself forvzard and down mtil it unites 
with the forward slope of the breaking wave, resulting in 
the formation of a vortex with a circulation around a 
cavity. VJhen the tip of the "jet” touches the forv?ard slope 
of the breaking wave, the resulting force forms a second 
weaker vortex in front of the first one. In this manner, a 
series of vortices yields a series of large, single pulse 
eddies, which contribute to the rapid decay of the wave 
during breaking. It has been found that the magnitudes of 
breaker vortices are generally correlated with a gradation 
from plunging to spilling breakers. 

Thornton et al. (1976) have observed a sloping profile 
toward the trailing edge of the wave. Secondary waves are 
frequently noted on the trailing edge both on the analog 
strip charts, Figure 1, and by visual observations in the 
field. The secondary waves have been observed to grow during 
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Figure 1. 



Strip Chart Record of Sea Surface Elevations 
and Particle Velocities. 
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the shoaling process, especially just prior to breaking 
when the primary waves are steepening significantly. The 
secondary V7aves are a result of the transfer of energy from 
the primary wave frequency to harmonics, and help maintain 
the potential energy across the surf zone. 

Gallagher (1971) has photographed secondary wave forma- 
tion occurring across reefs. He suggests that the secondary 
waves are brought about by a response to vertical accelera- 
tion near the front of a broken or breaking swell. When the 
amplitude and surface curvature of a shallow water wave ex- 
ceed certain values, the non-linearities of the vertical 
accelerations bring soliton wave forms into being. 

Figure 1 shows that the maximum particle velocity occurs 
near the time of maximum surface elevation. It can also be 
noted in Figure 1 that there are secondary velocity maxima 
which correlate with the secondary surface elevation maxima. 
The secondary velocity peaks also are a result of the non- 
linearity associated with the breaking wave. 

B. EIIERGY DENSITY SPECTRA OF BREAKING WAVES 

In order to better understand the PDFs as determined in 

this study, the energy density spectra measured by Galvin 

(1975) and Olsen (1977) were reviewed. The energy density 

spectrum is proportional to potential wave energy within a 

frequency band. The energy density spectrum, G (f) , indicates 

2 

how the variance, a , is distributed over frequencies, f, 

X 

where 
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00 



a 



2 

X 




G^(f)df 



( 1 ) 



(Kinsman, 1965), 

The energy density spectra, as estimated from the sea 
surface elevation records, were used to find the frequency 
of peak density and the occurrence of harmonics . The power 
spectra that were measured by Galvin showed a peak density 
over the frequency range of 0.061 to 0.092 Hz (16.4 to 10.9 
sec). First, second and third harmonics can be observed 
as shown in the spectra of surface elevation and horizontal 
velocities for 4 March 19 75. Additionally, subharmonics 
were observed at 0.011 and 0.922 Hz in the velocity spectra. 
Galvin (1975) has suggested that the observed low frequency 
sub-harmonic peaks may be due to edge waves in the surf zone. 

Thornton et al. (1976) note that the appearance of strong 
harmonics could be real due to secondary waves or artificial 
due to Fourier analysis of a peaked wave form. It is noted 
that secondary maxima of n are apparent in the strip chart 
record shown in Figure 1. It is felt that the secondary 
spectral peaks are real since secondary waves appear on the 
strip charts. 

G. PROBABILITY DENSITY FUNCTIONS OF PERIODS AlTD HEIGHTS 

Studies of the PDFs of period, height, and wave lengths 
have been conducted over the last 20 years. The majority 
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SPECTRAL DENSITY M /SEC ( M -SEC) 




Figure 2. Energy Density Spectra of n and v, 4 March 1975. 
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of these studies have been mainly concerned with theoretical 
or deep water aspects of the problem. 

The theoretical distribution of wave amplitudes for a 
narrow-banded Gaussian surface elevation in deep water, was 
shown by Longuet-Higgins (1952) to be a Rayleigh distribu- 
tion. Tayfm (1977), in studying the transformation of deep 
water waves to shallow water waves, showed that the Rayleigh 
distribution for wave amplitude was generally applicable to 
all bandwidths. Using the assumption that the wave height 
is twice the wave amplitude, the wave height PDF is then re- 
presented by: 



P(H)dH = ^ exp 



■H“ 

8a 



dH 



( 2 ) 



where a is the variance of the surface elevation. Using 
pressure records in the Gulf of Mexico, Longuet-Higgins 
(1975) observed that the Rayleigh distribution fit the ob- 
served distribution reasonably well. He found that there 
is a slight excess of waves with heights near the middle 
of the range and a deficit at the two extremes . Since much 
of the high frequency portions of the wave records are fil- 
tered out by the pressure transducer, Longuet-Higgins (1975) 
suggests that the narrow-band approximation may not be as 
applicable for the unfiltered records. In shallow x>?ater 
with much steeper waves, the distribution can again be 
expected to be less applicable due to the non-linearities 
which become increasingly important. 
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Since the statistics of H and T for a broad-banded 
spectrum vjere still somewhat unclear, Goda (1970) simulated 
wave profiles using an infinite series where the amplitudes 
were specified by various theoretical spectra and the phase 
was random with a computer. He then examined the simulated 
wave records for surface elevation, crest- to- trough wave 
heights and zero-up-cross wave heights. He found that for 
the zero-up-crossing determination of wave heights, the Ray- 
leigh distribution is a good approximation irrespective of 
the spectral width. 

Goda (1970) suggests that the deviation of the wave 
height distribution from the Rayleigh type in shallow v/ater 
should be considered separately from the change in the func- 
tional shape of the wave spectrum. It must be a highly non- 
linear phenomenon, which requires full analysis with the 
theory of non-linear wave spectra. 

Bretschneider (1959) investigated the PDFs of wave 
heights, H, wave lengths, L, and wave periods, T. He con- 
cluded from his empirical study that in deep water wave 
systems, the PDFs of wave heights and lengths can be approxi- 
mated by a Rayleigh distribution. Using the relationship, 
wave length proportional to period squared for deep water 
linear waves, he concluded that the PDF of periods squared 
is a Rayleigh PDF given by 



t2 

p (T) = 2.7-^ exp 


-0.675 


’I 








f 













where f is the mean period. 
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Using simulated wave data, Goda (1970) found that the 
crest-to-cres t periods squared were generally not well ap- 
proximated by the Rayleigh PDF. The mean and standard devia- 
tions of the periods by this method were found to be a 
function of the spectral width. He found that the periods 
determined by the zero-up-cross method are only in fair 
agreement with Brets chneider (1959). If the waves are of 
the type which have sharply peaked spectra, double peaked 
spectra, or a flat spectrum, Bretschneider ' s (1959) distri- 
bution was totally in disagreement with that of the simu- 
lated waves . 

Contrary to Bretschneider (1959) and Goda (1970), 
Longuet-Higgins (1975) noted that in general the deep water 
wave periods were bell shaped (but not Gaussian) . He defines 
the period PDF when calculated by the zero-up-cross method 
as , 

p(n) = ^^- 5 — j ; n = (t-<t>)/v<t> (4) 

2(l+n'^) ^ 

where <x> is the mean period and v is proportional to 
the spectral width where 

V = (y2/Po)^(<''^>/27r) (5) 

and and y 2 are the zero and second moments of the 

energy spectrum. 

Koele and de Bruyn (1964), Goda (1967), and Siefert 
(1970) have observed that the distribution of wave heights 
in shallow water does not correlate well with the Rayleigh 
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distribution. Kuo and Kuo (1974) have suggested that this 
is due to: 1) the non-linear effects of wave interactions 
yielding more larger waves, 2) the effect of bottom friction 
yielding a reduction in the low frequency components, and 
3) the effects of wave breaking, which would trmcate the 
distribution and transfer some of the kinetic energy to the 
high frequency components. Therefore the PDFs of wave 
heights with a certain intensity, may be considered to ap- 
proximate a conditional Rayleigh distribution truncated by 
the local breaking height. 

A zero-up- cross analysis was conducted on waves measured 
on a reef by Black (1978), who then compared the observed 
distribution with the Rayleigh, truncated Rayleigh, and 
the Weibull distributions. He observed that the truncated 
Rayleigh distribution fit well for those waves in the breaker 
zone when the heights were depth limited. The truncated 
Rayleigh distribution did not fit well for those waves 
measured either inshore or offshore of the breaker zone. 

He has found that a Ueibull distribution more closely approx- 
imated these measured distributions. 

It would be expected that the Weibull distribution would 
yield the better fit. The Rayleigh distribution is a func- 
tion of the variance of the data, whereas the Weibull dis- 
tribution is a function of higher moments about the mean. 
Neither theoretical distribution adequately describes ob- 
served distributions which are multi-modal. 

Black's (1973) observations of wave periods showed a 
distribution which was skewed to the right. This does not 
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fit well with the symmetrical Longuet-Higgins (1975) distri- 
bution. It was found that the period distribution more 
closely followed the Rayleigh distribution than the period 
square as proposed by Brets chneider (1959). It was found 
that the Weibull distribution most closely fit the empiri- 
cal. PDF, but there is considerable variation in the peaked- 
ness parameter. His attempts in relating the period statis- 
tics to the position on the reef were tinsuccessful. He 
feels that this is due to the non-linearities of the waves 
in shallow water. 

In this thesis, the empirical data collected for breaking 
waves were compared with the Rayleigh distribution. 

D. JOINT PROBABILITY DENSITY FUNCTIONS 

A study of the joint PDF for wave heights, n , and lengths, 
X, which can be viewed as heights and periods) was accom- 
plished by Brets chneider (1959). The general form of the 
joint PDF for two related variables (n,X) is the product of 

the marginal PDF (p(n)) and the conditional PDF (p(X)). 

n 

p(n,X) = p(n)*Pj^(X) (6) 

then 

CO 00 

p(n) = J p(n,X)dX = p(n) j p^(X)dX (7) 

o o 

Assuming that the marginal PDFs are Rayleigh, Bret- 
schneider (1959) examined the extreme cases of 0 and +1 cor- 
relation. For the case of zero correlation, the marginal 
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Rayleigh PDFs are given by 



p (n) 


= ^ n exp(-TTn^/4) ; 


n = H/H 


(8) 


P (A) 


= j A exp (-itA^/4) ; 


A = L/L 


(9) 



For zero correlation, i.e. total independence (Figure 3), 
the joint PDF is, 

p (n. A) = p (n) *p (A) (10) 

so that the joint PDF of heights and lengths is given as, 

2 

p(n,A) = ^ n exp (-Trn^/4) • A exp (-ttA^/4) (11) 

In terms of heights and periods, equation 15 is, 

p(n,T^) = 1.35 TT exp (-7rri^/4) • exp(-0.675 x”^) (12) 

where x = T/T. 

For the case of total dependence, the correlation coef- 
ficient is equal to one. For a correlation coefficient of 
one, all data points on a plot of joint Rayleigh PDFs fall 
on a 45 degree straight line passing through the origin, 
n=A=0. Assuming that a joint Rayleigh distribution is appli- 
cable, the following relation applies: 

p(n,A) = ^ n exp(-TTn^/4) • J X exp(-irA^/4) (13) 

which is plotted in Figure 3. 
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Figure 3. Statistically Independent Joint Rayleigh PDF. 
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Collins (1970) studied the PDFs of breaking waves de- 
rived from deep water waves having a joint Rayleigh distri- 
bution of heights, lengths and angles. He assumed that 
H = yd » where II is the breaker height, d is the breaker 
depth, and y is the breaking criterion assumed. The limit- 
ing factor of waves breaking in the surf zone yields a dis- 
tribution better approximated by a truncated Rayleigh dis- 
tribution. As expected, he fomd that the longest deep water 
waves were associated with the largest breaking waves. 

Collins (1970) concluded from the study of the joint 
PDFs of wave heights and lengths, that the effect of deep 
water angle of approach is relatively minor for breaking 
waves. The breaking wave steepness was found to be very 
sensitive to the choice of breaking criteria. 

A theoretical expression for the joint distribution of 
wave periods and amplitudes was derived by Longue t- Higgins 
(1975) for a narrow frequency spectrum 

2 

p(€,n) = — — n exp (-C^(l+n^)/2) (14) 

(2tt)'= 



where C and n ane the normalized amplitudes and periods 
as defined by. 



C = a/u_ 



% 



(15) 



and 



n = (x-<T>)/\;<T> 



(15) 



where \> is defined in equation 5. The wave amplitude 



22 



th 

is a , is the zero’^ moment of the energy spectrum, 

and <T> is the mean wave period. Longuet-Higgins (1975) 
compared this theoretical distribution of ocean waves to 
that of the deep water data obtained by Brets chneider (1959) 
and found good agreement. 

Chakrabarti and Cooley (1977) compared the theoretical 
joint distribution of Longuet-Higgins (1975) V7ith wave data 
recorded during a 1961 storm in the North Atlantic. The 
spectrum of these waves did not fall in the category of 
narrow-band. Through the comparison, it was found that 
there was a considerable agreement at the higher wave heights, 
but a definite trend away from the theoretical curve at the 
lower wave heights . 



23 



III. EXPERIMENT 



A. EXPERIMENTAL SITES 

Three different California experimental sites were used 
to obtain data on the three types of breaking waves--plung- 
ing, spilling, and collapsing. During the period 4 March to 
10 March 1975, plunging and spilling waves were measured by 
Galvin (1975) on Del Monte Beach, Monterey, California, 

This locale was chosen because the waves are highly refracted 
and directionally filtered, yielding narrowbanded swell off- 
shore which impinges almost perpendicularly on the beach. 

On 29 May 19 75 Galvin (19 75) conducted a second experi- 
ment on the Carmel River Beach, which is located approximate- 
ly five miles south of Monterey, California. This area is 
in an embaynent , where the waves are again highly refracted 
and directionally filtered. The beach at this site is very 
steep resulting in plunging and often collapsing and surging 
breakers. Due to the steepness of the beach and the rapidity 
of breaking of the shoaling waves, reflected waves are found 
to be present. 

Spilling breakers on a gently sloping beach were measured 
in March 1977 at Torrey Pines Beach, a site just north of 
La Jolla, California (Olsen, 1977). 

On 3 and 9 March 1978 data were collected by this author 
on Del Monte Beach, Monterey, California. On the first day 
of the experiment, very narrowbanded swell was observed 
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breaking on the shore, mainly of the plunging type. A storm 
arrived late on 8 March and continued through 9 March. V7inds 
of approximately 15 knots yielded a very confused sea. 
Breakers were of both the spilling and plunging type and 
were observed to break at various small angles relative to 
the beach. Sea surface elevations and velocities were col- 
lected for approximately three hours on each of these days 
on both the flowing and ebbing sides of high tide. 

B. INSTRUMENTS 

Capacitance type gauges were used to measure the change 
in the sea surface elevation as explained by Thornton et al. 
(1976). The wave gauges were constructed of 3/8 inch poly- 
propylene covered stainless steel rod. Calibration was 
accomplished in the laboratory statically by lowering the 
wave gauge into a deep water tank. Depth of immersion was 
then plotted against volts. 

The velocity of the water particles was measured with 
Marsh-McBimey Model 721 and 722 Electromagnetic Current 
meters. These flowmeters work on Faraday's principle of 
electromagnetic induction and are capable of recording the 
particle velocities in orthogonal directions in a plane. 

The meters were dynamically calibrated in the laboratory. 
Calibration was accomplished by oscillating the current 
meters in a water tank with a prescribed motion simulating 
the flow conditions which would be expected to be found in 
the surf zone. The flowmeter voltage output was found to 
be linear in the velocity ranges found in the breakers. 
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All instruments were placed on towers in the surf zone 
at low tide at the expected point of breaking of the waves 
for the subsequent high tide. Each tower held a capacitance 
wave gauge and an electromagnetic flowmeter. The flowmeter 
was placed below the wave gauge, approximately 0.5 m from 
the bottom, to insure constant submergence. Data were col- 
lected within an hour or two on both the ebb and flow sides 
of high tide. Analog recordings were made on both a strip 
chart and magnetic tape. 

C. DATA ANALYSIS 

All signals were digitized from their analog form by 
means of a Vidar Corporation 32-channel digital data acqui- 
sition system. The digitized data were then processed on 
an IBM 360/67 computer. Px.ecord lengths of approximately 
30 minutes from each data set were analyzed. Specific de- 
tails of the analysis can be found in Appendix A. 

The data were edited by excluding any value greater 
than six standard deviations from the mean and replacing 
that value with the previous value. The values were then 
linearly detrended to exclude the effects of the rising or 
falling tides. The mean and variance of sea surface eleva- 
tions were calculated for each record as well as the mean 
and variance of the heights, periods and maximum velocities 
all of which are defined on page 27. 

It was desired to only examine the sea-swell band of 
frequencies lying between 0.03-1.0 Hz (30 sec-1 sec). Hence, 
the data were passed through both a high and low pass digital 
filter. 
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The high pass filter had a cut-off frequency of 0.03 Hz 
(30 sec). The high pass filter used a Fast Fourier Trans- 
form (FFT) algorithm to obtain the Fourier spectrum of the 
entire 30 minute record. The Fourier coefficients correspond- 
ing to 0-0.03 Kz were then used to synthesize a low frequency 
time series which was subtracted from the wave record. 

The low pass filter was a 25 weight inverse transform 
filter (Davidson, 1970), as shovm in Figure 4. A cut-off 
frequency of 1.2 Hz vjas used with a terminal frequency of 
1.8 Hz. 

Wave heights were determined from the surface elevation 
record by means of the zero-up-cross method. The zero-up- 
cross technique defines a wave period as the interval between 
adjacent upcrosses. The maximiim and minimum of the surface 
elevation within the interval defines the crest and trough 
of a wave. Maximum onshore particle velocities were calcu- 
lated from the flowmeter records. Maximum onshore velocity 
was defined as the maximTom velocity amplitude occurring be- 
tween successive zero upcrossings as defined by the sea 
surface elevation record. Joint PDFs of wave height and 
period, wave height and maximum velocity, and maximum veloc- 
ity and period were calculated from the analyses. 

Heights and maximum velocities were normalized by divid- 
ing by their respective standard deviations. The periods 
were normalized by dividing by their mean period. The nota- 
tion used is 

H = H-/a^ 

V = V-/a^ 

T = T' /T 
27 
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Figure 4. Low Pass Digital Filter. 
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where heights (H' ) , maximum velocities (V), and periods 
(T') are the values as determined by the zero-up -cross 
method. T is the mean period, and and are the 

standard deviations of the surface elevations and particle 
velocities respectively. Histograms were calculated for 
the normalized heights, periods, and maximum velocities. 

The histograms were transformed into PDFs by dividing the 
number of value over an interval by the total number of 
values in the array and by dividing by the interval width. 
The mathematical expression for the exact PDF is 

TO “ 

p(x) = lim „ .lim ^ 

^ ^ Ax ^ 

where N is the number of times the variable occupies an 
interval, II is the total nxjmber in the record observed. 

Ax is the interval width and p (x) is the PDF. A restric- 
tion on the probability density is that 



(18) 



j p (x) dx = 1 . 
o 

Joint probability density functions of periods versus 
heights, velocities versus heights and periods versus 
velocities were calculated in a similar fashion. A 15 by 15 
grid was used to determine the frequency of values which 
simultaneously fill each joint interval. Mathematically 
the probability that the arrays x and y simultaneously 
assume values within the range (x, x+Ax) , (y, y+Ay) may be 
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obtained by taking the ratio of the total number of times of 
each occurrence within a bin, and the total number of 

observations, N . The exact probability density is approached 
as the number of observations approaches infinity the interval 
widths approach zero. The joint PDF can be defined as: 



p(x,y) = lim 



1 rlim 
X oo AxAy ^ N 

y oo 



N 

_iEZ 

II 



(19) 



with the restriction 

-f-oo 



II 



o o 



p (x,y) dx dy = 1 



The energy density spectrum (power spectrum) was obtained 
by means of a Fast Fourier Transform (FFT) algorithm. The 
energy density spectrum tells how the variance is distributed 
with respect to frequency. 

The spectral width parameter was calculated using the 
equation : 



£ 







2 



( 20 ) 



where mQ, m 2 , and m^ are the spectral moments (Cartwright 
and Longuet-Higgins , 1956). For e=0 the waves have only one 
frequency and the energy spectrum would be represented by a 
spike. As e-»-l the sea surface can be described as a broad- 
band process for which the distribution would be Gaussian 
and the spectrum would approach a horizontal line. Swell 
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can therefore be described as a narrow-band process and seas 
as broad-band. Due to the appearance of secondary waves, a 
narrow-banded deepwater spectrum of waves broadens in the 



surf region. 

The correlation coefficients of the joint values were 
calculated to obtain a measure of the amount of linear rela- 
tionship of random pairs of values of periods and heights , 
velocities and heights, and periods and velocities, 



where x and y are the sample means, and y^ are 

the sample standard deviations. 
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IV. DISCUSSION 



Miller (1976) concluded that breaker shapes form a con- 
tinuous gradation of shapes from spilling to plunging. Each 
breaker type undergoes the same processes, but differs only 
in the magnitude of the vortices generated. The probability 
density fmctions appear to show a similar gradation. Two 
simplified cases, for spilling and plmging, were chosen to 
represent extreme examples and are used for all discussions. 
The 4 March 1975 data most simply represent the plunging 
case and the 16 March 1977 data represent the spilling case. 
The results for all data are given in Appendix C. 

In both examples for spilling and plunging breakers, the 
offshore waves appeared as narrow band swell. A spectral 
analysis of the surface elevation records in the surf zone 
yielded spectral width parameters ranging from 0.94-0.97 
(see Appendix A), suggesting a broad band spectrum in the 
surf zone. The broad spectral width can be attributed to 
the non-linear wave-wave interactions yielding harmonics in 
the breaking zone. 

In order to determine the origin of the large spectral 
width parameter, c was calculated in the ranges of 0-0.25 
Hz, 0-0.50 Hz, 0-1.0 Hz and 0-nyquist frequency. The analy- 
sis revealed that large values of the spectral width para- 
meter occurred in all cases. In the case of 0-0.25 Hz v;hich 
contains the first three harmonics, e was greater than 0.7, 
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which was the case with the smallest e . It is concluded 
that the change from a narrow deep water spectrum to a 
broad-band shoaling spectrum is primarily due to generation 
of secondary waves at harmonic frequencies. 

A. PERIOD PROBABILITY DENSITY FUNCTIONS 

Bretschneider (1959) made the assumption that the wave 
lengths were Pv.ayleigh distributed; using this assumption 
and the linear wave theory relationship 

L-^t2 ^ (22) 

the wave periods squared would be Rayleigh distributed in 
deep water. Brets chneider' s analysis can be extended to 
the shallow water case using the shallow water approxima- 
tion for wave length 

L = /gK T (23) 

Assuming the wave lengths are Rayleigh distributed, the 
periods in shallow water would be expected to be Rayleigh 
distributed. 

The period PDFs for plunging and spilling examples are 
shown in Figures 5 and 6 respectively. The Rayleigh PDF has 
been superimposed on the period PDFs for comparison. The 
formulas used for the Pvayleigh PDF for various parameter 
transformations are given in Appendix A. In both the spill- 
ing and plunging wave examples, it is noted that the period 
PDFs only loosel^^- follow the Rayleigh PDF. The Rayleigh PDF 
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Figure 5. PDFs of Period and Height, 16 March 1977. 
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Figure 6. PDFs of Period and Height, 4 March 1975. 



overpredicts the wave periods greater than the- Rayleigh mode 
value, and underpredicts the observed periods less than the 
peak Rayleigh value. The tr-uncation of the long period 
waves is probably due to the larger period waves breaking 
farther offshore and then reforming to be recorded by the 
stationary wave gauge as higher frequency waves. 

The spilling example for 16 March 1977, shown in Figure 5 
has the mode lying close to the maximum of the superimposed 
Rayleigh PDF. The peak of the breaking wave distribution 
is located at the first harmonic of the offshore wave period, 
which is attributed to the secondary waves. A peak corre- 
sponding to the peak offshore wave period is absent. 

The plunging waves, as represented by the 4 March 1975 
record in Figure 6, are multimodal. Since the Gaussian, Gram 
Charlier, Rayleigh and Weibull PDFs are unimodal, none of the 
theoretical PDFs will adequately describe multi-modal PDFs 
for this particular plunging surf condition. 

Two strong peaks are foimd to occur in the plunging 
breaker period PDF. The longer period is found to correlate 
well with peak period of the wave spectrum. The shorter wave 
period appears at approximately the first harmonic. Since 
the center period of the bin is used to represent the inter- 
val, the magnitude of the peaks and the proximity to har- 
monics as represented by the labels is a function of the 
interval width and the interval cutoffs. A 30 bin PDF was 
calculated for better resolution to determine if the peaks 
were harmonics and Figure 7 shows that they are. 
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Figuie 7. PDFs of Period and Height, 30 bin, 4 March 1975, 



The waves periods as determined by a zero-up- cross of 
the surface elevation record for spilling conditions yields 
a poor indication of the peak of the energy spectrum which 
corresponds to the mean offshore wave period. In the spill- 
ing breaker case, the mean period corresponded to the first 
harmonic of peak in the energy spectrum. For the plunging 
breakers, the mean period corresponded to a valley in the 
period PDF between the primary and first harmonic periods. 

The appearance of strong peaks at one half the offshore wave 
period is indicative of the importance of non-linearities 
in the form of secondary waves in the surf zone. 

B. HEIGHT PROBABILITY DEIISITY FUNCTIONS 

Wave height PDFs for the plunging and spilling wave ex- 
amples are shown plotted against the theoretical Rayleigh 
PDF in Figures 5 and 6. The height PDF for spilling waves 
is not as easily described as the period PDF. The height 
PDF is loosely described by the Rayleigh PDF. The major 
height peak is centered at 2.2. The Rayleigh PDF overesti- 
mates the density of the small wave observations and tends 
to underestimate the density of the large waves. The heights 
are truncated at approximately six standard deviations. Since 
the wave gauge is stationary and therefore the depth of the 
water is essentially fixed, this truncation occurs due to 
the breaking of the large waves seaward of the wave gauge. 

The calculated mean height value on 16 March 1977 lies at 
2.7, which actually corresponds to a dip in the PDF relative 
to the peak intervals on either side. 
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The mean height observed in the 4 March 19 75 plunging 
breaker PDF was calculated to be at 2.4. The observed PDF 
is fomd to be strongly bimodal with peaks at 1.8 and 3.4. 

The strong bimodality is attributed to the very narrow band 
incident wave system and the generation of secondary waves 
at the first harmonic. 

C. JOINT PROBABILITY DENSITY FUNCTIONS (T vs H) 

The joint PDF of periods and heights for the plunging 
and spilling examples are contoured in Figures 5 and 6; the 
distributions for all cases are found in Appendix C. 

It is noted that the longest wave periods are associated 
with the greatest wave heights suggesting strong correlation 
between periods and heights . The correlation coefficient 
was calculated for all cases and varied between 0.60 and 0-.80. 
Bretschneider (1959) examined the two special cases of joint 
Rayleigh distributions with correlations of 0 and +1. Per- 
fect correlation between two variables would result in 
correlation coefficient of 1.0 and a straight line in the 
joint PDF. 

The joint Rayleigh PDF for the case of a correlation co- 
efficient of zero is plotted in Figure 3. Zero correlation, 
or statistical independence, results in the joint PDF being 
spread out as compared with the straight line for perfect 
correlation. The amount of spread is a measure of the reduc- 
tion in correlation. 

Figure 5 is an example of p(H,T) for plunging breakers 
and shows tv;o regions of high density. The bimodality of 
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p(H) and p(T) are reflected in the p(H,T). The maximum peak 
period, which corresponds to the peak in the wave spectriom, 
is highly correlated with the primary peak of wave heights. 

The second high density region corresponds to the first 
harmonic of the primary wave period. The joint PDF suggests 
strong wave-wave interactions in a plmging wave system re- 
sulting in energy transfer to harmonics. 

Two distinct contour groups are no longer apparent in 
Figure 6 for the spilling case, but rather one long contour 
is found. Since the height PDF and the period PDF are only 
weakly bimodal, it follows that a single contour group 
should be found. The area of maximum joint density is foimd 
to be associated with the first harmonic of the offshore 
wave period, and the mode of the observed breaker heights. 

Greater density was fomd at short periods and heights 
for all records than predicted by the Rayleigh PDF. A high 
frequency wave riding on a long period wave causes an addi- 
tional zero crossing as the longer wave passes through zero. 
Most previous comparisons of heights or periods of waves have 
been done using pressure sensors which greatly filter the 
high frequencies due to hydrodynamic attenuation. Capaci- 
tance wave staffs and electromagnetic flowmeters have good 
response times and thus measure high frequencies. Therefore, 
even though the data were filtered at 1 IIz , all records showed 
peaks at the shortest periods and heights. 
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D. VELOCITY PROBABILITY DENSITY FUNCTIONS 



The velocity PDFs were calculated as the maximum onshore 
velocities within the period segments determined by the zero- 
up-crosses of the surface elevation record as shown in Fig- 
ures 8 and 9. Using this technique can result in both posi- 
tive and negative maximas for onshore particle velocities 
(i.e. , d^v/dt < 0). 

The velocity PDFs for spilling breakers and plunging 
breakers are shown in Figures 8 and 9. The theoretical 
Rayleigh PDF is superimposed on the same graphs for compari- 
son. It is apparent that neither the spilling or plunging 
case is adequately described by the Rayleigh PDF. The veloc- 
ity PDF for spilling breakers appears to be symmetrical, with 
the mode value falling at the bin with center velocity of 
1.5. There also appears to be a small peak value at 0.3. 

When compared with the plunging case, it is found that the 
peaks occur at the same normalized velocities. The plunging 
wave PDF is bimodal and has relatively few occurrences of 
negative onshore velocities. 

E. JOINT PROBABILITY DENSITY FUNCTIONS (V vs H) 

The joint probability distributions of velocities and 
heights are shown in Figures 8 and 9. The results show high 
correlation between velocities and heights. The high posi- 
tive correlation suggests that the larger maximum onshore 
velocities are correlated with the larger wave heights. 

Figure 8 shows the plunging breaker case. Two peak areas 
are found to occur in the joint PDF which reflects the 
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Figure 8. PDFs of Velocity and Height, 16 March 1977. 
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Figure 9. PDFs of Velocity and Height, 4 March 1975. 



bimodality of the individual PDFs. The joint PDF for the 
plxjnging waves suggests a highly correlated case (0.68) in 
which the largest onshore velocities occur due to waves of 
greatest height. 

The spilling case as shown in Figure 9 also shows high 
correlation (0.70). For this example, there is only one 
major peak which corresponds to the individual PDFs which 
are essentially unimodal. The peak area of the joint PDF 

A 

has been shifted toward higher waves and greater velocities 
than the plunging case. 

F. JOINT PROBABILITY DENSITY FUNCTIONS (T vs V) 

Figures 10 and 11 show joint PDFs of periods and veloc- 
ities in which the periods were determined from the zero-up - 
crosses of the wave gauge record and the coincident maximum 
onshore velocities are determined from the flowmeter record. 
The results are similar to the joint PDF of period and 
heights. Spilling waves, as contoured in Figure 11, show one 
main contour group. The maxim\im spilling probability is 
found at the intersection of period 0.84 and velocity 1.5. 

The plunging waves as found in Figure 10 reveal two contour 
peaks. One large contour group is associated with the point 
representing the offshore period 1.56 and a velocity of 1.5. 

The second large contour is a function of the first harmonic 
wave period of 0.84 and a velocity of 0.3. 

The negative values of the maximum velocity indicate 
that, during a wave period defined by the zero-up-crosses 
of the surface elevation, the maximum velocity can be negative. 
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Figure 10. PDFs of Period and Velocity, 16 March 1977. 
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Figure 11. PDFs of Period and Velocity, 4 March 1975. 



The joint PDF of wave height and velocity shows that the 
negative velocities are associated with the low wave heights. 
The joint PDF of velocity and period shows that small nega- 
tive velocities are correlated with the short periods. 
Therefore, it is concluded the negative velocity maxima 
are associated with short periods and low wave heights. The 
secondary waves occur after the passage of the primary wave 
crest and have a peak near the trough of the primary wave; 
the velocities associated with the primary wave would be 
negative maximum in this region so that the measured veloc- 
ity, which is a sum of the primary and secondary waves, can 
be small or negative. 
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V. CONCLUSIONS 






Waves in the surf zone are a highly non-linear process 
which is evidenced by the appearance of secondary waves. 

The secondary waves appeared in the PDFs at the first har- 
monic of the primary period. 

The joint probability density fxmctions of the breaking 
wave properties show high correlations. The correlation 
factor ranged between 0.60-0.80 for periods and heights 
and suggests that the greater the wave period, the higher 
the breaker height will be. High correlation between the 
periods and velocities and the heights and velocities sug- 
gests that the maximum onshore velocities are a fmction 
of wave periods and breaker heights. 

Breaking waves may be viewed as a gradation between 
two extremes, spilling, and plunging breakers. The differ- 
ence between the spilling and plunging breakers becomes 
apparent in the joint PDFs. The spilling breakers yield a 
joint PDF which has only one high density contour region, 
whereas the plunging breakers show two high density contour 
groups. The high density contour groups appear to be most 
associated with the period PDF. 

The joint PDF (p(T,H)) of spilling breakers has a maxi- 
mum density at a wave period equal to one half the peak 
period of the wave energy density spectrum, i.e. the first 
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harmonic. The corresponding velocities (p(v,T)) or heights 
(p(H,T)) were likewise fomd to be mimodal with a peak den- 
sity which is driven by the first harmonic wave period. 

Plunging breakers, at the opposite end of the gradation, 
yield a bimodal joint PDF. The two major contour groups 
appear to be a result of the bimodal density distribution 
of the periods. The two peak densities of the period PDF 
are representative of the peak of the wave density spectrum 
and a first harmonic. Due to the high correlation of periods 
and heights, the corresponding PDF of wave heights appears 
to be bimodal. The joint PDF of period and velocity and the 
joint PDF of velocity and height suggests high correlation 
between these variables. The maximum onshore velocities 
appear to be correlated with both wave periods and wave 
heights . 

A trmcated P^ayleigh PDF was found to loosely approxi- 
mate the empirical period and the height PDFs of spilling 
breakers. This truncation seems to be the result of long 
period waves with their associated larger breaker heights 
breaking prior to the arrival at the wave gauge. The mode 
of the spilling wave period distributions corresponds to the 
first harmonic of the peak of the energy density spectrum. 

The smoothed appearance of the spilling wave periods PDF, 
relative to the plunging case, appears due to the smearing 
which occurs with the appearance of harmonics greater than 
the first. The velocities PDFs of spilling waves, where the 
velocity maximum is determined as in Section D, yields a 
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symmetrical distribution with some negative onshore veloc- 
ities present. Since the Rayleigh distribution is undefined 
for negative values, the resulting fit of observed data with 
the theoretical Rayleigh distribution is poor. 

The PDFs of the plunging breakers are poorly described 
by either the Rayleigh PDF or other existing theoretical 
PDFs such as the Weibull PDF. Theoretical PDFs are unimodal, 
whereas the empirical PDFs under plunging conditions are 
strongly bimodal. The high correlation found in the joint 
PDFs for wave heights and maximum velocities also results 
in a bimodal distribution both in the velocities and the 
heights PDFs . 

The mean wave period as obtained by the zero-up-cross 
method in the siirf zone is a poor indicator of the mean 
offshore wave period. Spectral analysis of the breaking 
waves shows an energy peak corresponding to the primary off- 
shore wave period, but the period PDFs of breaking waves 
have a strong peak at the first harmonic of the wave spectra 
peak. Plunging waves had high density peaks in the period 
PDF corresponding both to the peak of the wave spectrxam and 
to its first harmonic. The result is that the mean of the 
breaking waves falls between the primary and the first har- 
monic wave period. The mean period of breaking waves there- 
fore is a poor estimation of the mean period of the incident 
waves . 
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APPENDIX A 



STATISTICAL PARAMETERS CALCULATED 



1. Due to the variation of sampling frequencies between the 
different experiments, it follows that some of the other 
statistical parameters will also vary. A summary of these 
parameters follows where f^^ is the Nyquist frequency, Af 
is the resolution. At is the interval between samples and 
T is the length of sampled record in minutes. 



DATS 




%(Hz) 


Af (Hz) At(s 


amp le / s e c ) 


T (min) 


4 


MARCH 


1975 


2.5 


6.1x10"'^ 


0.200 


27.3 


5 


MARCH 


19 75 


2.5 


6.1x10"'^ 


0.200 


27.3 


6 


MARCH 


19 75 


2.5 


6.1x10'^ 


0.200 


27. 3 


8 


MARCH 


19 75 


2.0 


4. SxlO"'^ 


0.256 


35.0 


29 


MAY 


19 75 


2.0 


4.9x10'^ 


0.256 


35.0 


16 


MARCH 


1977 


2.0 


4.9x10'^ 


0.250 


34. 1 


17 


MARCH 


19 77 


2.0 


4.9x10"^ 


0.250 


34.1 


19 


MARCH 


1977 


2.0 


4.9x10"^ 


0.250 


34. 1 


8 


MARCH 


19 78 


2.0 


4.9x10"^' 


0.250 


34.1 


9 


MARCH 


1978 


2.0 


4.9x10"^ 


0.250 


34. 1 


2. 


Tables of 


calculated statistical 


. parameters 


follow 




n = 


sea : 


surface elevation (m) 








V = 


instantaneous 


particle velocities (m/sec) 



VM = maximum onshore particle velocities (m/sec) 

LR = length of record for calculated statistics (min) 
e = spectral width parameter 
H = wave height (m) 

T = wave period (sec) 
r = correlation coefficient 
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APPENDIX B 



MATHEMATICAL DERIVATION OF NOPxMALIZED RAYLEIGH DISTRIBUTIONS 

Wave heights, periods and velocities have been compared 
with the Rayleigh distribution 

A a2 

p(A) = — 7T exp (- -^) 

4a 2a^ 



where A is the non-normalized value. Since all parameters 
are normalized in this paper, the following conversion was 
made : 



let A' = A/a 
then 



,dA 



p(A') = p(A) (^, ) A=A’a 
p(A') = p(A) = p(A)a 



for wave heights 



p(H) = exp(-%(^)2) 






exp(-%(^)2) 

n n n 



for maximrun velocities 



P(V) = -^2 exp(-%(^)2) 

40v V 



/V X V / 1 / V X 2x 

p(— )= TTT- exp(-%(j^) ) 






V 



for periods 



P(T) = 



4a 



^ exp(-%(^)2) 
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where 



a = f/2 



then 



P(r) 




exp(-%(^)2) 



P(“) = r exp(-%( 
T T 
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Hi iH 



APPENDIX C 



PROBABILITY DENSITY FUNCTIONS 

The following tables consist of both the normalized 
joint probability density functions, and the individual 
probability density functions as labeled. The PDF tables 
coincide with Figures 5 thru 11, respectively. 

The top table is the joint PDF. The top line and the 
left column represent the center period, height, or 
velocity point of each of the 15 bins. The numbers in 
the other columns and lines represent the frequency of 
occurrence in the respective bins. 

The bottom table represents the individual PDFs as 
labeled. The top line of the individual PDFs is the 
center point of the individual bins. The second line 
represents the frequency of occurrence in the bin. 
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PDFs of Period and Height, 4 March 1975. 
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PDFs of Period and Height, 16 March 1977. 
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